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Abstract 

In this paper we propose an approach to prepare GHZ states of an arbitrary 
multi-particle system in terms of Grover's fast quantum searching algorithm. 
This approach can be regarded as an extension of the Grover's algorithm to 
find one or more items in an unsorted database. 
PACS number (s): 32.80.Pj, 42.50.vk 
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The original motivation to prepare three-particle entangled states was from the obser- 
vation by Greenberger, Home, and Zeilinger (GHZ) that entanglement of more than two 
particles leads to a conflict with local realism for nonstatistical predictions of quantum me- 
chanics [1]. This is in contrast to the case of experiments with two entangled particles testing 
Bell's inequalities, where the conflict only arises for statistical predictions [2]. However, up 
to now it has become a well known fact that the significance of GHZ states is far beyond the 
original motivation. In fact, a GHZ state plays a key role in quantum teleportation [3] which 
lies in the heart of quantum information. Recently, GHZ states of three particles [4] have 
been produced experimentally. On the other hand, the field of quantum computing [5-8] 
has undergone a rapid growth over the past few years. Simple quantum computations have 
been performed experimentally. Especially, two great quantum algorithms, Shor's factoring 
algorithm [9] and Grover's quantum searching algorithm (GQSA)[10] have been experimen- 
tally realized for those cases of a few qubits [11,12]. So it becomes an interesting subject to 
investigate applications of these quantum algorithms. In this paper, we make use of Grover's 
quantum searching algorithm to prepare GHZ states of an arbitrary multi-particle system. 

We consider a system consisting of n qubits, it has N = 2 n possible independent states. 
We symbolize these states as | 0), | 1), • • • , | TV— 1). In Grover's quantum searching algorithm, 
there are two elementary unitary operations. One is the operation M performed on a single 
qubit that is represented by following matrix: 



One can perform the transformation M on each qubit independently in sequence to change 
the state of the system. This operation is called as Walsh-Hadamard transformation [13] 
which is expressed as 



If the initial state of system is in | 0), then the resultant state after performing transformation 
(2) is the identical superposition of all N states. This is a way of creating a superposition 
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with the same amplitude in all iV states. 

The other elementary operation in the Grover's algorithm is the selective rotation of the 
phase of the amplitude in certain states. The transformation describing this operation for a 
n-qubit system is of the following form: 



X = 











(3) 



where ip , (fii, ■ ■ -(fiN-i are the arbitrary real numbers. 

Making use of the above two elementary unitary operations, we can construct the fol- 
lowing iterative operation: 



q = -wxiw-^x* 



(4) 



where Xo andx 71 are two selective rotation transformations of the phase of amplitude, and 
given by 



Xo 
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We now use the Walsh-hadamard transformation (2) to initialize the system. After 
initializing, the state of the system becomes 



N-l 

*> = ^|°> = ^X>> 



(5) 
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Then, we successively perform the operation (4) on above initialized state. After j 
iterations, we obtain the following state, 

N-2 



Q 1 I ^o> 



kj-= (| 0)+ | N - 1)) + lj-= V 1 i), 



(6) 



with 



N 



kj = \/? sin [(2j + 1} 9] ' lj = V ATT2 C ° s[(2j + m 
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where 9 is defined by sin 2 9 = 2/N. 

If we select the iterative times j to satisfy: 
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then we have = 0, kj = yj y. From Eq.(6) we therefore obtain the following n-bit GHZ 
state: 
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At this moment, if we perform the following selective rotation of the phase of amplitude on 
the GHZ state (9), 
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we can get another n-bit GHZ state: 



^>6=^(|0>-|iV-l». 



(11) 



Unfortunately, that j determined by Eq. (9) is not always an integer. In what follows we 
will present two methods to avoid this problem. 

The first approach is to prepare a new initial state. We note that the above Grover's 
operations are related to the initialized state. We only slightly modify the initialized state 



determined by Eq.(5), then we can get satisfactory result. In order to see this, we choose 
ji = [j + 1], an integer that slightly larger than j. By using j 1 = — | and sin 2 #i = 2a\ 
to determine a\. Clearly a\ is slightly smaller than 1/y/N. If we can produce the following 
new initial state: 



ti) = ai (\0)+\N-l)) + 



\-2a\ 
N -2 
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then just by performing the former GQSA ji times on Eq.(ll), we immediately achieve the 
GHZ state in Eq.(9). The new initial state in Eq.(12) can be easily obtained through the 
following Grover operation of arbitrary phase [14]: 



Qi = -Wx^w-'x?, 



(13) 



where Xo an d X 13 are given by, respectively, 
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By runing the operation (13) on the initialized state (5), we prepare the following state: 



N-2 



Q l \^ ) = k(\0}+\N-l}) + lJ2 



(14) 
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where 



k = 2(e to -l)_+^ + (e fa -1)(N- 2) 
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Comparing the above state with the state (12), we can find that k — a\ and 

. ai Ny/N + N-A 
cosa = cos/3 = 2 (N-2) ^ ' 



Since a\ is slightly smaller than 1/v A, a and (3 are two small angles. Therefore, we realize 
the new initial state (12). 

The second method is to slow down the speed of the very last iteration. In this case we 
choose jo = [j], an integer being smaller than j, then run GQSA j times on the initialized 
state (5), it becomes: 

JV-2 



| Vo> = 0)+ I N-l)) + l j0 ^= £ | i), 

* * i=i 



(16) 



where 



k j0 = Yy sin [(2 jo + 1)9], l j0 = yj jj—^ cos[(2j + 1)9], (17) 
where sin 2 9 = 2/N. 

After above operation we again perform the following operation with an arbitrary phase 
on the state (16): 

Q 2 = -WxtW-'x* (18) 

then we get 

N-2 

Q 2 Q jo I V>o> = fc'(| o)+ | iv - l)) + r \i) (19) 



1=1 



where fc' and I' are given by 



k' = 2 ^-^ +jV e^„ + ^- 1 ^- 2 ' i ,-„ (20) 



= 2 <^>e^„ + (e» -!)("- 2) + AT 
Av^ J Av^ J 



From Eqs.(18)-(21) we see that the GHZ state (9) can be realized, if we properly modulate 
and (p to satisfy /' = in Eq.(19). This condition I' — can be satisfied by requiring that 
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(22) 



N — A N-2kl 




(23) 



cosy? = 



Finally, we simply analyze two special cases of three qubits and two qubits. For the case 
of three qubits, n = 3 and N = 8. we only need to perform one time GQSA opteration on 
the initialized state (5), then we can immediately reach the common GHZ state: (| 000}+ | 



Grover operation with 7r/2-phase, i.e., the transformation in Eq.(13) with a = [3 = ir/2, on 
initialized state in Eq.(5), then we can get the widely useful Bell's states. 

In summary based on the GQSA we have proposed an approach to produce GHZ states 
of arbitrary number qubits. This method can be also used to exactly search one or more 
items from a database. Grover's original quantum searching algorithm is actually a quantum 
amplitude amplification . It can amplify the probability of the searched items in a database 
very approaching to one, but not certainly equal one, depending on the probability of the 
searched items. However, Grover operations of arbitrary phase can successively amplify or 
decrease a small probability of the searched items. We have combined above two algorithms 
to realize a successful and certain quantum searching to a database. In this sense, we can say 
that the method suggested in present paper is the extension of Grover's quantum searching 
algorithm. 
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For the case of two qubits, n 



2 and N — 4. We may perform one time 
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